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We  consider  the  asymptotic  distribution  of  functionals  associated  with 
normed  sums  of  random  plane  convex  sets.  Methods  involving  symmetric  statistics 
and  weak  convergence  of  stochastic  processes  are  used  to  examine  area  and 
perimeter  in  particular. 
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EXPLANATION 

Minkowski  addition  and  scalar  multiplication  of  subsets  of  JR^  are  defined 
respectively  by 

\ + K2  = {kx  + k2  : kx  e KL,  k2  e 
aK1  = {ak^  : k^  e K^ } . 

If  Xl,X2,...  are  random  subsets  of  TR^  then  it  is  of  interest  to  study  the 
behavior  of  their  averages 

X = - [X  + •••  + X ] . 
n n 1 n 

Following  on  earlier  work,  we  consider  the  behavior  of  area  and  perimeter  of 

X when  the  X.  lie  in  the  plane, 
n 1 r 

Potential  applications  for  this  work  include  the  modelling  of  biological 
growth  phenomena  and  techniques  for  two-dimensional  image  processing. 
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ASYMPTOTIC  AREA  AND  PERIMETER  OF  SUMS  OF  RANDOM  PLANE  CONVEX  SETS 

Richard  A.  Vitale 

1 . Introduction. 

If  Xj/Xj,...  are  random  subsets  of  ]R^,  then  we  may  form  their  successive 

averages,  X = — (X  + •••  + X ],  under  Minkowski  addition  and  inquire  into  their 
n n 1 n 

asymptotic  behavior.  A first  step  in  this  direction  was  taken  in  [1]  where  a strong 
law  of  large  numbers  for  sets  was  derived.  Our  purpose  here  is  to  look  at  some 
distributional  considerations  associated  with  this  convergence  in  the  case 
when  the  X^  are  valued  in  the  compact,  convex  subsets  of  the  plane.  Certain  of  our 
results  evidently  have  extensions  to  higher  dimensions  but  we  have  chosen  the  restriction 
d = 2 for  uniformity  of  exposition.  The  assumption  of  convexity,  on  the  other  hand, 
is  central  to  most  of  the  discussion. 

The  paper  is  divided  as  follows.  In  Section  2 we  set  down  some  notation  and 
preliminaries  including  a statement  of  the  strong  law.  Section  3 presents  a central 
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2.  Preliminaries. 


Let  K denote  the  collection  of  non-empty,  compact,  and  convex  subsets  of  TR  . 

With  the  introduction  of  the  Hausdorff  metric 

d(Kl(K2)  = inf {e  > 0 : ^ £ K2  + cB,  K,,  < + eB^  < 

K may  be  regarded  as  a separable,  locally  compact  metric  space.  Here  B is  the  closed 
unit  disc  and  scalar  multiplication  and  Minkowski  addition  are  defined  as  usual : 

eK  = {ek  : k e K} 

K1  + K2  = + k2  : kx  ( K1#  k2  e . 


The  norm  of  K is  defined  by 


= maxi  k : k e K} 


or,  equivalently,  d ( { 0 } , K) . 

To  each  K,  we  assign  a support  function  given  by 

(2.1)  s(0,K)  = max{k-e(0)  : k f K)  e(0)  = (cos  0,  sin  0) 

0 t [0, 2wJ  . 

The  map  K -*  s(-,K)  uniquely  embeds  K in  the  space  C[0,2n).  Linear  structure  as 
well  as  distance  are  preserved: 

(2.2)  s(-,aK)  = as(-,K)  a 0 

(2.3)  s ( • , K + K2)  = s ( • ,K  ) + s(-,K2) 

(2.4)  . d(Kl(K2)  = ||  s ( - ,KX>  -s(-,K2)|| 

= max{  | s (0  ,K  ) - s(0,K2|  : 0 e [0, 2tt  ] > 

IIkII  = || s ( - ,kj  || 

We  shall  regard  a random  set  X as  map,  measurable  in  the  Borel  sense,  from  an 
abstract  probability  space  into  A'.  In  order  to  formulate  the  strong  law  of  large 
numbers,  let  us  first  define  the  expectation  of  a random  set.  By  virtue  of  the  inequality 

| s (0  ,X)  | < ||x|| 

Es(0,X)  is  well-defined  if  e||x||  < ■». 

Definition.  Let  e||x||  < ■».  Then  the  expectation  of  X,  written  EX,  is  that  element 
of  K which  uniquely  satisfies 
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3 . Weak  Convergence  of  the  SupjKirt  Function  Process 

To  each  random  set  X is  associated  its  support  function  s(-,X),  which  may  be 
regarded  as  a random  element  of  Cl0,2n).  Using  (2.2)  and  (2.3)  we  have 


, 

;(0,x  ) = ~ l s (0 , x . ) 


i=l 

and  the  difference  process 

1 n 

(3.1)  A (6)  = s(0,X  ) - s(0,EX.)  = - J (s  (0,X.  ) - ES(0,X.))  . 

n n in.,  l r 

1=1 

Theorem  1 . Let  X^.X^....  be  iid  random  sets  with  e||xJ|2<  °°.  Then,  as  n -*■  “, 

47  An(0)  converges  weakly  to  the  C[0, 2t] -valued  Gaussian  process  6(0)  with 

EA (0)  = 0 

Cov(A(0l)  ,A(02>)  = Cov(s(0  ,X.),s(02,Xi))  . 

Proof.  The  existence  of  EX^  is  certainly  ensured  by  the  condition  in  e||xJ|2<  <». 

In  addition,  this  requirement  together  with  the  representation  (3.3)  of  6^(0)  as  a 

sum  of  iid  random  values  implies  that  the  finite  dimensional  distributions  of  47  A (0) 

n 

converge  to  the  appropriate  Gaussian  limits. 

It  remains  to  show  tightness  of  the  measures  on  C[0,2n]  associated  with  the 
47  A (•).  Accordingly,  we  verify  the  bound 


(3.2) 

(12,  p.  95)).  We  have 


E I 47  6 (0  ) - 47  A (6, ) 1 < e||x.||  (0,  - 0,)' 

n 2 nl  — " l"  2 1 


E[47  A (0  ) - 47  6 (0,))2  = E [s  (0  , X . ) - s (0  , EX . ) - s(0,,X.)  + s(0,,EX.))2 
n 2 nl  2 i 2i  li  li 

< E [s  (0  , X . ) - s (0  ,X.))2  . 

— 2 i 1 l 

Now  an  easy  application  of  the  Cauchy-Schwarz  inequality  together  with  (2.1)  provides 
the  general  bound 

|s(02,K)  - s(01(K'  j < ||k||  |02  - 01| 

which  yields  (3.2). 


As  an  example,  let  us  consider  the  homogeneous  line  segment  model  which  is  constructed 


as  follows.  Let  a be  uniformly  distributed  on  (0,2ti)  and  let  L,  independent  of  a, 

2 

be  a non-negative  random  variable  with  EL  < Then  X is  taken  to  be  the  line 

segment  of  length  L centered  at  the  origin  and  inclined  at  an  angle  a to  the 

horizontal.  It  follows  that  s(0,X)  = ^ |cos(0  - a) | and  EX  is  the  disc  of  radius 
EL 

— centered  at  the  origin.  The  covariance  structure  of  s(0,X)  and  hence  of  A is 

TT 

given  by 

£i=M.(i.MHh|j.(a]2. 

Theorem  1 and  (2.4)  imply  a rate  estimate  for  the  convergence  Xn  -*■  EX. 

Corollary . With  the  assumptions  of  the  theorem,  i^n  dfX^.EX)  converges  in  distribution 
to  ||  A ( - ) ||  - 

Unfortunately  at  this  stage  we  know  of  no  methods  for  investigating  ||A(-)||  in 
even  the  simplest,  non-trivial  cases  such  as  the  example  just  given. 

The  situation  for  bounded  linear  functionals  is  direct.  Let  T be  a map 
taking  K into  IR1  such  that  Tfa^  + eK2)  = aTO^)  + BT(K2>  a, 8 >_  0 and 

|t(K)  | C ||  K ||  for  some  C > 0 and  all  K f K. 

Corollary ■ With  the  assumptions  of  the  theorem,  a bounded  linear  functional  obeys  a 
central  limit  theorem:  Jn  (T(X^)  - T(EX)]  converges  in  distribution  to  an  N(0,Var  T(X^)) 
variable. 

The  perimeter  functional  provides  an  example  for  this  convergence.  Here 
2ti 

per(K)  = / s(0,K)d0  so  that  /n(per(X  ) - per(EX))  converges  in  distribution  to  a normal 

0 " 2i  2n 

variable  with  mean  0 and  variance  Var  per (X . ) = / / Cov(s (0 ,X. ) ,s (y,X. ) )d9dy  . Xn 

1 0 0 1 1 

the  homogeneous  line  segment  model  we  can  evaluate  this  expression  directly  by  noting  that 
per(XJ  H 2L^  (the  factor  2 enters  since  a line  segment  is  regarded  as  a degenerate 
polygon)  and  so  Var  perfX^)  = 4 Var  L^ • 

Related  linear  functionals  are  T (K)  = s(0o,K),  the  extent  of  K in  the  direction 
0Q,  and  T(K)  = s (Oq, K)  + s(00  + n,K),  the  width  of  K in  the  direction  0g. 


Cov(s(0  + h,X),s(0,X))  = 


EL 
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4.  Area. 


Let  us  first  remark  that  the  continuity  of  the  area  functional  ensures  that 

A(X  ) -*  A(EX)  a.s.  under  the  conditions  of  the  strong  law  result.  Evidently  the  mode 

n 

of  this  convergence  must  be  investigated  by  means  other  than  those  in  the  preceding 
section  since  area  is  not  linear  in  the  Minkowski  sum. 

Instead  we  use  the  representation 


(4.1) 


A (X  ) 
n 


1_ 

2 . 

n 3 


l A (X  . ; X ) . 
,k  3 K 


Here  A(X.;X.)  = A(X.)  while  A(X.;X  ) (j  ^ k implicitly)  is  the  symmetric  mixed  area 
3 1 3 3 k 

of  X.  and  X (14)).  We  shall  assume  throughout  that  E ||  X . || 4 < " which  ensures  a 
3 K 3 

bounded  second  moment  for  each  term  of  (4.1)  (note  0 A(K  ;K  ) ^h||kJ|  Hk^H  ). 

The  form  of  (4.1)  suggests  the  use  of  standard  techniques  for  studying  symmetric 
(or  U-)  statistics  (as  discussed,  for  instance,  in  von  Mises  (8),  Hoeffding  [7j, 
Filippova  (5J , Rubin  and  Sethuruman  (91).  With  some  labor  in  the  adaptation  this  can 
be  done.  However,  we  shall  take  a different  approach  which  is  based  on  orthogonal 
expansions  (a  general  discussion  of  this  method  will  appear  in  [10]). 

Since 

EA  (X  ) = --  [ (n2  - n)  EA  (X  . ; X ) + nEA(X.)l 
n n2  3kg 


we  have 

A (EX  ) = lim  EA  (X  ) = EA  (X  ; X ) 

1 n 12 

so  that 

(4.2)  6 = A(X  ) - A (EX, ) 

n n 1 

= A'  l (A  (X  ; X ) - EA  (X  . ; X ) I 
n j^k  3 3 K 

+ ~ l (A  (X  . ) - EA(X  )) 
n 3 3 3 

+ - (EA (X, ) - A (EX, ) ) . 

n l i 
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We  define 

°2  = Ev  1Ev  A<X  :X.  ) - EA(X.;X.  )]2,  j * k 
a . A "1  K 1 K 

3 k 

2 

which  is  finite.  The  convergence  properties  of  6 will  depend  on  whether  o is  zero 

n 

or  strictly  positive. 

Theorem  2.  Under  the  assumptions  made  above,  6 can  display  two  types  of  convergence: 
n 

2 1/2  2 

(i)  If  a >0,  then  n 6 converges  in  law  to  a mean  0,  variance  4o 

n 

normal  variable. 

2 

(ii)  If  a =0,  then  n6^  converges  in  law  to  a variable  of  the  form 


l c (Z*  - 1]  + EA (X  ) - A (EX  ) 
V=1 


where  the  Z^  are  independent  standard  normal  variables  and  the  c^  form  a square 
summable  sequence  described  below. 

Proof.  Before  considering  the  two  cases  separately,  we  make  the  following  observations. 
Under  either  normalization  the  contribution  of  the  second  term  in  (4.2)  is  asymptotically 
negligible.  Moreover  the  bias  effect  of  the  third  terra  persists  only  in  the  second 
case.  Accordingly,  we  focus  on  the  first  term 


(4.3) 


S = l I A (X  . ;X  ) - EA(X.,X  ))  . 

n n j*k  3 k 3 k 


Each  term  of  this  expression  is  symmetric  in  its  arguments  and  has  bounded  second  moment. 
This  is  sufficient  to  conclude  a denumerable  expansion  (convergent  in  mean  square)  of 
the  form 

A (X  , ; X ) - EA  (X  . ; X ) = Y c V>  (X.)*  (X.  ) 
jk  jk  ^vvjvk 

where  (X^)  = 6^  ((3,  p.  1087)).  We  then  have 


(4.4) 

(4.5) 


o = e[a(x.:x  ) - ea (x . ; x ) l = Y c y2 
lk  3 k *■  v v 

V 

2 r 2-2 
o = ) C 

u V V 
V 


where  <p ^ *=  EV v (X^ ) - Finally  we  have 
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Sn  = ~T  II  7v(Xj)#v,\> 
n 3*k  v J 

= K l/v'V^v^1  • 

n v ytk  J 
r -2 

Case  (i) . With  the  insertion  of  the  null  term  L c V ((4.4))  we  have 


nl/2sn  = n‘V2  l %(<  I ‘W  - V Wv<\>  - %> 
v k 


(V\,(X.)  - V ) W (X,  ) - Vj] 


+ 2 (n  - l)*’ 


I ‘W  ' *vj} 


l CV>  I <VXj>  - V(VV  ' V1 

v j^k 


2 (n  - 1 )n  T c 7 (^  (X.)  - <p  ) . 

v v ^ v i v 
V j 


Using  the  orthonormality  of  the  , one  can  verify  directly  that  the  variance  of  the 

v 

first  term  tends  to  zero  with  increasing  n (it  is  in  fact  bounded  above  by 
4n  2 ) £ • T*le  secor>d  term  is  essentially 


V = 2n’1/2  lev  l {V  (X.)  - V ) 


u V V v VI  V 

V 3 J 


for  which  we  consider  a finite  truncation 


V = 2n  l c V l V X.)  - V . 

nN  v v t v j v 

v=l  3 J 


Again  using  the  properties  of  the  <p  , one  can  find 

EV  = EV  „ = 0 
n nN 


E(V  - V J2  < 4 l CV 
" ""  ~ v=N+l  V V 


N 2_2  N 

EV2  = 41  7 C V - ( Y c V>  ))  . 
nN  V v *■,  v v 

v=l  v=l 


Hence  V -*■  V in  mean  square  uniformly  in  n.  Since  the  inner  sum  of  (4.6) 
nN  n 

is  over  iid  elements  with  zero  mean  and  finite  variance,  we  conclude  that,  as  n ■*  *»,  V 
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converges  in  law  to  a zero  mean  normal  variable  with  variance  given  by  (4.8). 

2 

With  increasing  N,  (4.8)  converges  to  4o  . This  is  sufficient  to  conclude  the  result 

(see,  for  instance,  (2,  Theorem  4.2)). 

2 

Case  (it)  . Here  o = 0 so  chat  by  (4.5)  = 0 for  all  v.  Again  we  use  a trunca- 


tion argument.  Setting 


V = nS  = - l c l V>  (X.)l?  (X,  ) 
n n " v V j^k  v > v k 


V = - l c l * (X.)*>  (X.  ) 
nN  n , v .i,  v 1 v k 
v=l  J*k  J 


EV  = EV  = 0 
n nN 


oo 

(4.9)  E [V  - V 1 2 < 4(  "In-2  j c2 

n nN  “ l2'  V=N+1  V 

As  before  (4.9)  provides  a uniform  bound  on  the  truncation  error  (in  mean  square).  Now 


we  may  re-write 


S'-'vl’-i  j*v] 

N i 2 N 

= 1 c | -J-r  l f (X.)  - J c M V»2(X.)  . 

v=l  v 1 nV2  i v 3 ' v=l  V " j V J 

N 

By  the  law  of  large  numbers,  the  second  term  converges  to  £ c . Setting 

-1/2  r V=1  V 

V =n  ) f (X.),  we  have 

n\l 


.VI 

1 


EY  — 0 
nv 

EY  Y =6 
nv  nvi  vp 


so  that  (Y  Y ) converges  in  law  to  (Z_,...,Z  ) where  the  Z are  independent 

nl  .nN  IN  v 

N(0,1)  variables.  Hence,  as  n -*■  V converges  in  law  to  a variable  of  the  form 

nN 

N o 

l c (Z2  - 11 

L.  V V 
V=1 
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which,  with  increasing 


N,  converges  to 


l c \Z  - 1) 

v=l  v v 


and  we  are  done . 

As  an  example,  let  us  again  use  the  homogeneous  line  segment  model.  Here 


A (X  . ;X  ) = | sin  (a  . - a )j  . 

j k 2 j k 1 


We  have 


EA (X . ;X  ) = 

1 k 


(ELk>2  2 = (KV2 

2 tt  7i 


and 


so  that 


EEL  L.EL 

A (X  ; X ) = -V1  !=  -V5- 


Xk"'-3  k 


2 it 


2 (EL^  V 7 
o = __  var  I- 

n 


nV2  A(X  )■ 
n 


If  Var  I.  > 0,  then  a >0  and  the  central  limit  theorem  holds.  That  is, 
k 2 


(Eh.) 


converges  in  law  to  a normal  variable  with  zero  mean  and  variance  4 
In  the  other  case,  let  us  suppose  for  concreteness  that  L 1 . Then  we  have  the 
ordinary  Fourier  orthogonal  expansion 


A (X  . ; X ) - EA  ( X . ; X ) = — I sin  (a  . - a,  ) - E sin  (ci . - a ) ] 

j k 3 k 4 ' 3 k 1 3k 


l 


i 


v=2.  ,4  , . . . 1 - v 


cos  va  . cos  va,  sin  va  . sin  va, 
i k . i y 

77  • 77 — + ~7T^  ■ 7T  ~ 


Since  EA(XJ  = 0,  case  <ii)  provides  the  convergence  in  law 

- v 1 - — S '<  - 1 4 - 11  - ; - 1 - v - , J 

V— 2 , 4 , . . . 1 - V v-  2,4,. 


where  we  have  used 


v=2,4,...  v 


— ([6,  #6.1.37])  and  indicated  by 


2 - 2 

{W  = Z + Z } a collection  of  independent,  mean  2 exponential  variables. 
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